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Semi-leptonic .D+(1968) decays as a scalar meson probe 
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The unusual multiplet structures associated with the light spin zero mesons have recently attracted 
a good deal of theoretical attention. Here we discuss some aspects associated with the possibility 
of getting new experimental information on this topic from semi-leptonic decays of heavy charged 
mesons into an isosinglet scalar or pseudoscalar plus leptons. 
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, I. INTRODUCTION 

<\ 

j- , ■ For the first years after the quark model was accepted it was believed that the lightest scalar meson should be a 

' quark-antiquark composite with mass value similar to those of the tensor and axial vector mesons. In particular, an 
occasionally discussed light, broad "sigma meson" was not expected to exist. However, more recent work has provided 
evidence for such a particle as well as for other similar light scalars. (Some characteristic references are [l[ - (38[.) 

In fact there seem to be enough scalar candidates to fill up two different nonets. A model including these states, 
with also two nonets of low l yin g pseudoscalars in order to form chiral multiplets has been studied in some detail; 
plj . [22| . [27| . [30| . 32|-[36j], [38| . Note that chiral symmetry is the exact symmetry of QCD with massless quarks. 
Adding "soft" quark mass terms results in "sigma models" which give many reasonable low energy predictions. In 
i ' the models just mentioned one chiral (i.e. containing both scalars and pseudoscalars) nonet is supposed to represent 
states with a quark-antiquark substructure while the other nonet is supposed to represent states with a two quark - 
two antiquark substructure. The physical states are suitable linear combinations. 

On the experimental side of the subject, information on the light scalars has often been extracted from study of 
pion pion and other scattering processes. Another way is to search for scalar resonances explicitly in particle decay 
ly-j processes. Recently, the CLEO collaboration has reported [39| good evidence for the scalar / (980) in the semi-leptonic 
(y-) ■ decay of the D+(1968) meson. Since there is more phase space available, it may be possible to find other scalar iso- 
singlet states in this and similar semi-leptonic decays of heavy mesons. There are also isosinglet pseudoscalar states 
like the r] and 7/(980) which can be studied and in fact have been already reported in the decays of the D+(1968). 

As a possibly helpful adjunct to future work in this direction we will, in the present paper, make some theoret- 
ical estimates of the semi-leptonic decay widths of the D+(1968) into the four scalar isosinglet states and the four 
pseudoscalar isosinglet states which are predicted in the chiral model mentioned above. 

In section II we discuss the hadronic "weak currents" which are needed for the calculation. These are mathematically 
given by the so-called Noether currents of the sigma model Lagrangian being employed. We work in the approximation 
where renormalization of these currents from the symmetry limit are neglected. This means that there are no arbitrary 
parameters available to us. Nevertheless there are some subtleties. To explain these we build up the model in three 
stages rather than just writing the final result immediately. 

In section III we give a detailed description of the calculation of the partial decay widths from the currents discussed 
in section II. For this purpose we also use information on the scalar and pseudoscalar meson masses and mixings 
obtained in [3rj| . 

A short summary and discussion is given in section IV. 

In Appendix A we briefly discuss the well known K£3 decay which has the same general structure as the semi- 
leptonic D s decays. 
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II. HADRONIC CURRENTS IN VARIOUS LINEAR SIGMA MODELS 

These models give the usual " current algebra" results near the threshold of pion-pion scattering but also yield some 
additional interesting features away from threshold. 

A. Chiral SU(3) model 

The usual chiral nonet M(x) realizing the qq structure of the pseudoscalar and scalar mesons is schematically 
written with chiral SU(3) indices displayed as: 

Ml = ( (?M )t 74 i±^ 9aA; (1) 

where a and A are respectively flavor and color indices. For clarity, on the left hand side the undotted index 
transforms under the left SU(3) while the dotted index transforms under the right SU(3). The decomposition in 
terms of scalar and pseudoscalar fields is: 

M = S + i<t>. (2) 
Using matrix notation (e.g. M b — > M b ) the Noether vector and axial currents read (see for example Appendix A 

of El), 

V„ = i^d^ + iS d„ S, 

A M = Sd^-^d^S, (3) 

The axial symmetry breaking is measured by the vacuum value of S: 

S = S+<S>, <S b a >= a a 5 b a , (4) 

where the normalization is u\ + 012 = F^ ss 130.4 MeV and a± + 03 = Fk ~ 156.1 MeV. Note that the overall 
normalization constant for gives the correct value for the ordinary electromagnetic current. This determines the 
normalization for the weak currents in the SU(3)l x SU(3)# symmetry limit. For the vector currents this amounts 
to an implementation of the "conserved vector current hypothesis" introduced for beta decay many years ago [4lj |. 
Such an approximation is well known not to be as good for the axial current case, but may at least furnish an order 
of magnitude estimate. In detail, with the usual SU(3) tensor indices, the currents read: 

Ka = X € + iS c a C S b c + i{a a - a b )d^l 

Ka = S% S b c + (a a + a b )d^ b a , (5) 

For example, the relevant hadronic current needed to describe the semi-leptonic decay K + — > 7r° + e + + v e is 

VZ^^d^l+i^-a^Sl (6) 
A relevant application of this formula is given in Appendix A. 

B. SU(3) M - M' model 

For this model we first introduce another chiral field, constructed out of two quarks and two anti-quarks as: 

M^ b = e acd e b ^(M^(M^. (7) 
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This object has the form of a "molecule" made of two M's. Alternatively one can schematically make two quark- 
two antiquark states denoted by and from a diquark combined with an anti-diquark in two different ways 
[2l| . One might as well consider the most general linear combination of M^ 2 \ and as a field representing 
an object, M' made from two quarks plus two antiquarks. M' has the decomposition, 



M' = S' + : 



(8) 



Then the Noether currents involve the sum of pieces constructed from the unprimed fields and from the primed 
fields. The latter take the form, 



s' a c a M ^ - C ^ s' c b + {p a + (3 b )d^' a b , 



(9) 



wherein, 



S' = S'+ < S' >, 



< s' b a >= pjl 



(10) 



The total currents are denoted as: 



V'jtotal) 
A b m (total) 



rb 

[ b 

ft a 



V, 



fj,a ' 



A'! 



(11) 



In contrast to the chiral SU(3) model above, all the primed and corresponding unprimed fields mix to give physical 
fields of definite mass. As a simple example, the transformation between the physical tt + and 7r' + fields and the 
original fields (say </> + and </>' + ) is (32| : 



r'+ 



cos t) 7: 
sin ft. 



— sin v m 
cos On 



fx 



(12) 



which also defines the transformation matrix, R n . 

The pion decay constant as well as (formally) the decay constant for the much heavier 7r(1300) particle are defined 
by the part of the axial current linear in the fields: 



Al^total) = F w a M 7r+ +F„ l d fi ir'+ + •••, 

= ( a i + ol 2 ) cos 6^ -(/3i + j3 2 ) sin ft„, 
F„> = (ai +a 2 ) sin 6„ + (A+^ 2 ) cos 6 n . 



(13) 



The angle 9^ depends on the detailed dynamics. |32| 

In what follows it will be useful for us to specify the mixing matrix for the four isoscalar scalar mesons in this 
model. A basis for these states is given in terms of the four component vector / = (f a , fa, f c , fj) where, 



fa 
h 
fc 
fd 



si + sj 

V2 



ss. 



S[ + S 2 



V2 



S3 nnnn. 



(14) 



In the above, the quark content is indicated on the right for convenience. Note that s stands for a strange quark 
while n stands for a non-strange quark. However these basis states are not mass eigenstates. Again, the detailed 
dynamics of the model is required to specify this. For typical values of the model's input parameters (see [36l ) the 
mass eigenstates make up a four vector, F = Lq 1 f with, 
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(O = 



0.601 0.199 0.600 0.489 

-0.107 0.189 0.643 -0.735 

0.790 -0.050 -0.391 -0.470 

0.062 -0.960 0.272 -0.019 



(15) 



The physical states are identified, with nominal mass values, as 



F = 



/o(600) 
/o(980) 
/o(1370) 
/o(1800) 



(16) 



It will also be interesting for us to give the typical result of the model for the mixing of the four isoscalar pseu- 
doscalars. The analogous basis states are: 



'la 



v/2 



rjb = 4>l ss, 

, i 1 +^ 2 
r?c = f= — nsns, 

A'3 



Vd = 4>£ nnnn. (17) 

For typical values of the model's input parameters (see [36]) the mass eigenstates make up a four component vector, 
P = i?o 1 r ! with, 



P : 



7?(547) 

ri(958) 
77(1295) 
77(1760) 



(18) 



(These identifications correspond to the favored scenario discussed in section V of [3g]). The dynamically determined 
mixing matrix is then: 



-0.675 0.661 -0.205 0.255' 

0.722 0.512 -0.363 0.291 

-0.134 -0.546 -0.519 0.644 

0.073 0.051 0.746 0.660 



(19) 



C. Hybrid M -M' model with a heavy flavor 



As recently discussed in (38|, the case of three flavors is special in the sense that it is the only one in which a two 
quark-two antiquark field has the correct chiral transformation property to mix (in the chiral limit) with M. In order 
to respect this property when a heavy meson is included in the Lagrangian, we should demand that "heavy" spin 
zero mesons be made of just one quark and one antiquark. In a linear sigma model the kinetic term would then be 
written as: 



C = -^Tr\d,Md^) - irr^M'^Af't), (20) 

where the meaning of the superscript on the trace symbol is that the first term should be summed over the heavy 
quark index as well as the three light indices. This stands in contrast to the second term which is just summed over 
the three light quark indices pertaining to the two quark - two antiquark field M' . Since the Noether currents are 
sensitive only to these kinetic terms in the model, the vector and axial vector currents with flavor indices 1 through 
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3 in this model are just the same as in Eg. dTTI) above. However if either or both flavor indices take on the value 4 
(referring to the heavy flavor) the current will only have contributions from the field M. This should be clarified by 
the following example, 



ifa d^ a c + is c 4 s?, 

Std^^-fad^S*. (21) 

Here the unspecified indices can run from 1 to 4. This equation is correct by construction but does not tell the whole 
story since the connection between the fields above and the physical states involves, as in the preceding cases, the 
details of the non- derivative ("potential") terms of the effective Lagrangian. 



V^{total) = = 
AUtotal) = A% = 



III. DIFFERENT SEMI-LEPTONIC DECAY MODES OF THE D+(1968) 

The initial motivation for this work was the recent experimental discovery (39| of the semileptonic decay mode, 



#+(1968) -> /o(980)e + ^ e , (22) 
in which the /o(980) was identified from its two pion decay mode. 

A relevant generalization is to consider other scalar isosinglet candidates than the /o(980). For example the SU(3) 
M - M' model contains four different isoscalar scalars, Fj. In addition, there are four different isoscalar pseudoscalars 
in that model, P^. Here we shall calculate the predictions of that model for all eight of these decays in the simplest 
approximation. This should provide some useful orientation. In fact there are no parameters which have not already 
been determined in the previous treatment (3(| of the model. 

The usual weak interaction Lagrangian is, 



C = ^j={J- W+ + J+ W- ) , (23) 



wherein, 



J+ = ^7^(1+75)017 + ^(1 + 75)^- (24) 
Here the column vectors of the quark fields take the form: 



U 



u 
c 
t 



D = 



d 
s 
b 



(25) 



and the CKM matrix, V is explicitly, 





Vud 


Vus 






V = 


v cd 


Vcs 




(26) 




v td 


Vts 


v tb 





A picture describing the relevant D s decays is given in Fig. ([T]). 
The corresponding semi-leptonic decay amplitudes are thus, 



amp(D+(p)) 



{ Pi(q)\ +e+(k) + l/(l)) % G f { <P i (q)\V* i (total)\Df( J >)> \ 

xs(l)7^1 + 7 5 )«(k), (27) 
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s _ 

" R(q)orF(q) 



vJI) 



FIG. 1: D a decay. 

where the spinor u(k) represents the outgoing e + and represents the outgoing v e . The relevant hadronic 
operators can be rewritten in terms of the mass eigenstate scalar isosinglets and the pseudoscalar isosinglets using 
Eqs. HU) and $TE$i as: 

V^(total) = iD+Z </>! + ••• 

= iD+Y i (R ) v d l .P j +--- (28) 

3 

A^total) = I); i)„ Sf ■ ■■■ 

= n; £ /•;, • ••• (29) 

3 

The transposed matrices Lq and i?o are given in Eqs. (| 15[) and (|19[) respectively, based on a typical numerical 
solution for the model parameters (36[. Next the amplitudes are given by, 

amp(Df(p) -> | ^ | + e +(fe) + „ e (|)) = ^|V CS | j^J. } (p„ + + 7is)«(k), (30) 

The squared amplitudes, summed over the emitted lepton's spins, are then, 

G fI^| 2 ^2 { ((ijj^)a } t 2fc ■ (P + 9)1 ■ (p + ?) - * • *(P + ?) 2 ], (31) 

wherein m e has been set to zero except for the overall 1/m 2 factor. 
This yields the unintegrated decay width, 

dv _Gi\v cs nm^f\ m{Ds) \±_ (32) 



d|q| 12tH» \ ((L ) 2l ) 2 J qo 
For integrating this expression we need, 

m 2 (D s ) — mf 



2m{D s ) 

where rrii is the mass of the isosinglet meson Fi or Pj and also the indefinite integral formula, where x = |q| 



(33) 



— \Jx 2 + m? - ^m 2 xJ x 2 +mf + ^mfln(x + \J x 2 + m?). (34) 



^x 2 + m 2 4 

Table U summarizes the calculations of the predicted widths, for D+ decays into the four pseudoscalar singlet mesons 
{Vi = ?7(547), i] 2 — 77(982), 773 = r/(1225), 7/4 = 7/(1794). Notice that the listed masses, m, are the "predicted" ones in 
the present model) and leptons. 
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m, (MeV) 


(Bo)* 


(q ma x)i (MeV) 


(MeV) 


553 


0.661 


906.20 


4.14 x 10" 11 


982 


0.512 


739.00 


7.16 x 10~ 12 


1225 


-0.546 


602.74 


2.57 x 10~ 12 


1794 


0.051 


166.31 


2.65 x 10~ 17 


TABLE I: pseudoscalars. 


rrii (MeV) 


(Lo)'2i 


(? ranl ), (MeV) 


Ti (MeV) 


477 


0.199 


933.23 


4.56 x 10" 12 


1037 


0.189 


710.79 


7.80 x 10 -13 


1127 


-0.050 


661.30 


3.62 x 10" 14 


1735 


-0.960 


219.21 


3.85 x 10~ 14 



TABLE II: scalars. 

Tabic HIl with the same conventions, summarizes the calculations of the predicted widths for decays into the 
four scalar singlet mesons ((/i, /a, ■ • •) = (<r, /o(980), • • •)) and leptons. 
Experimental data exist for only three of these eight decay modes. 

T(L»+ -> r]e + v e ) = (3.5 ± 0.6) x 10~ n MeV 
T{D+ ^ri'e + v e ) = (1.29 ± 0.30) x 10" 11 MeV 

r(L>+ ->■ foe+v e ) = (2.6 ±0.4) x 10~ 12 MeV (35) 

It is encouraging that even though our calculation utilized the simplest model for the current and no arbitrary 
parameters were introduced, the prediction for the lightest hadronic mode, T(Dj — » rje + v e ) agrees with the measured 
value. In the case of the decay Df — >• rje + v e the predicted width is about 30% less than the measured value. For the 
mode Df — > / (980)e + ^ e our predicted value is about one third the measured value. Conceivably, considering the 
large predicted width into the very broad sigma state centered at 477 MeV, some of the higher mass sigma events 
might have been counted as /o(980) events, which would improve the agreement. It would be very interesting to 
obtain experimental information about the energy regions relevant to the other five predicted isosinglet modes. 

Furthermore, these width predictions are based on Eqs. (fT5|) and (fl"9]) corresponding to particular choices for the 
quark mass ratio A3/A1 and the precise mass of the very broad 11(1300) resonance. Varying these within the allowable 
ranges gives rise to the allowed range of predictions displayed in Figs. ([2]) and ([3]). One can see that raising m[LT(1300)] 
and/or lowering A3/A1 yields better agreement for the predicted semi-leptonic decay width of the /o(980). Clearly, 
the simple model here provides reasonable estimates for the semileptonic decay widths of the 15^(1968). 

IV. SUMMARY AND DISCUSSION 

We saw that the partial widths for semi-leptonic decays of the .D+(1968) into isoscalar scalar singlets and pseu- 
doscalar singlets plus leptons could be well estimated in a simple model where the hadronic current was taken to be 
the Noether current associated with a minimal linear sigma model. 

The agreement between experiment and theory was better for the decays into the rj and rf than for the decay into 
the /o(980). The former involve the hadronic vector current, which is "protected" according to the conserved vector 
current hypothesis, while the latter involves the "unprotected" axial vector current. 

Clearly it would be interesting to try this technique for other semi-leptonic decays of charmed mesons and also for 
bottom mesons. We considered the case when the charged lepton was e + rather than the cases of /i + or t + . In those 
two cases an additional form factor as in the calculation of the K£3 decay discussed in Appendix A should be taken 
into account. 

Information about the scalars, involving however more work for disentangling the effects of the strong interaction, 
can also be obtained from the non-leptonic decay modes of the charm and bottom mesons. A treatment of Dg — > 
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FIG. 2: Starting from the upper left and proceeding clockwise: The dependences of the pseudoscalar partial widths on the 
current quark mass ratio A3/A1 and on the value of the 11(1300) mass. 



/o(980) + 7r + has already been carried out [42j]. The study of the decay like B+ —> scalar + e + + v e might be useful 
for learning about mixing between a cc scalar and the lighter three flavor scalars. 

A straightforward, but not necessarily short, improvement of this calculation would be to include both vector and 
axial vector mesons in the starting Lagrangian from which the currents are calculated. 
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Appendix A: K£3 decay 

As an illustration of Eq. © we consider the matrix element, between an initial K + state with 4-momentum k and 
a final 7r° state with four momentum p, of the strangeness changing vector current Vjfi , 

< 7r°(p)\V*\K+(k) >~ f+(t)(k+p)^ + f-(t)(k-p) M , (Al) 
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FIG. 3: Starting from the upper left and proceeding clockwise: The dependences of the scalar partial widths on the current 
quark mass ratio Az/A\ and on the value of the 11(1300) mass. 

where t = -(k — p) 2 . 

The first term of Eq. © contributes at tree level to the /+ form factor while the second term contributes to the 
/_ form factor. These two contributions are illustrated in Figs. 0^, and|4]3 in which the W boson which is connected 
to the leptonic current acts at the points X. Here we are evaluating this matrix element in the framework of the plain 
SU(3) linear sigma model in which, furthermore, the vector and axial vector mesons have not been included. 



u°(p> K+M 



K + (k-p) 

(a) (b) 
FIG. 4: K£3 decay hadronic current. 



According to the usual Feynman rules, 



/+ 



1 

VI' 
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_ 1 oa-oi ml -m 2 
f ^ L o 3 + ai J[ ml-t J ' {A2} 

wherein m K denotes the mass of the strange scalar particle contained in this model. Furthermore, the explicit form 
of the K kit coupling constant in the model was used in the expression for /_ [2l| . Notice that the first bracket in the 
equation for /_ evaluates to about 0.16 and that the physical kappa mass is about 800 MeV in the plain SU(3) linear 
sigma model. 

It is interesting that this decay allows one to learn something about the properties of the kappa meson. For this 
purpose it is necessary to use the process where a final fj, + is observed rather than a final e + . That is because the 
contribution of f-(t) to the decay width is proportional to the final lepton mass. Of course the effect of the if* (892), 
which contributes importantly to the f+(t) form factor should also be included to get increased accuracy. 
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